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I.  INTRODUCTION 


The  ability  to  compute  the  base  region  flow  field  for  projectile  configur¬ 
ations  using  Navi er-Stokes  computational  techniques  has  been  developed  over 
the  past  few  years.  This  capability  is  most  important  for  determining  aerody¬ 
namic  coefficient  data  and  in  particular  the  total  aerodynamic  drag.  The 
total  drag  as  described  in  this  paper  consists  of  the  pressure  drag  (excluding 
the  base) ,  viscous  drag  and  base  drag  components.  At  transonic  and  low  super¬ 
sonic  speeds  the  base  draq  component  is  the  major  contributor  to  the  i„tal 
aerodynamic  drag. 

The  majority  of  base  flow  calculations  to  date  have  modeled  the  base 
region  as  a  flat  solid  surface.  Many  of  the  actual  configurations  have  some 
form  of  base  cavity.  General  opinion  has  been  that  the  inclusion  of  a  base 
cavity  or  modifications  to  the  interior  cavity  of  a  projectile  base  would  have 
little  or  no  effect  on  the  overall  flight  performance  parameters. 

The  M325  projectile  under  certain  conditions  is  expected  to  be  aeroballis- 
tically  similar  to  its  parent  configuration  the  M4R3A1.  The  M32S  has  an 
aluminum/steel  base  which  is  configured  as  a  flat  cavity  (standard).  A  recent 
product  improvement  program,  undertaken  to  reduce  the  production  costs  and 
improve  shell  integrity,  resulted  in  the  design  of  a  new  base  configuration 
made  from  steel  and  containing  a  dome  cavity.  A  series  of  aeroballistic 
tests1  were  conducted  in  the  Transonic  Range  Facility  of  the  US  Army  Ballistic 
Research  Laboratory  (3RL) ,  to  determine  any  difference  in  the  aeroball i sties 
which  may  occur  between  the  standard  and  dome  base  configurations.  As  a 
result  of  these  tests,  differences  in  aerodynamic  performance  were  found  to 
exist  between  the  two  rounds.  The  most  significant  changes  in  the  aerodynamic 
data  were  in  the  lift  and  static  moment  coefficients.  The  drag  was  found  to 
be  reduced  by  a  few  percent  with  the  dome  configuration  having  the  lower  drag. 

A  computational  study  was  undertaken  to  determine  the  ability  of  the 
present  Navi er-Stokes  codes  to  predict  these  differences  and  to  further  under¬ 
stand  the  fluid  dynamic  behavior  which  can  account  for  such  small  changes. 
The  use  of  Navi er-Stokes  codes  can  provide  a  detailed  description  of  the  flow 
field  associated  with  the  M825  configuration  as  well  as  the  integrated  aerody¬ 
namic  coefficients.  The  work  reported  here  has  been  accomplished  using  an 
axi symmetric  base  flow  code.  Numerical  computations  have  been  completed  for  a 
Mach  number  range  of  0.8  <  M  <  l.S.  The  computed  aerodynamic  drag  show  the 
same  effect  as  the  experimental  data,  that  is,  a  small  reduction  in  the  total 
aerodynamic  drag.  Qualitative  features  of  the  computed  flow  field  are 
presented  in  the  form  of  mach  contours  and  particle  traces. 


II.  GOVERNING  EQUATIONS  AMO  SOLUTION  TECHNIQUE 
1.  GOVERNING  EQUATIONS 

The  complete  set  of  time-dependent  thin-layer  Navi er-Stokes  equations  is 
solved  numerically  to  obtain  a  solution  to  this  problem.  The  numerical  tech¬ 
nique  used  is  an  implicit  finite  difference  scheme.  Although  time-dependent 
calculations  are  made,  the  transient  flow  is  not  of  primary  interest  at  the 
present  time.  The  steady  flow,  which  is  the  desired  result,  is  obtained  in  a 
time  asymptotic  fashion. 
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Tin1  t  imp-dependent ,  thin-layer,  Navi  er-Slok.es  equations  written  in  strong 
conservation  law  form  for  the  a »i  symmetri  c  formulation*-’  are: 


is .  «  1 6K  t  -  „  _i_ «  . 

6  a  6C  6  c  Re  6C 

The  general  coordinate  transformati ons  are  defined  as: 

f,  =  £(x,y,z,t)  is  the  longitudinal  coordinate 

n  =  n(x,y,z,t)  is  the  circumferential  coordinate 

r,  =  c(x,y,z,t)  is  the  near  normal  coordinate 

t  =  t  is  the  time. 


(1) 


The  vector  q  contains  the  dependent  variables  [p,  pu,  pv,  pw,  e]  and  the  flux 
vectors  E,  G  contain  terms  which  arise  from  the  conservation  of  mass,  momen¬ 
tum  and  energy  in  the  three  coordinate  directions.  The  source  vector  H,  in 
Equation  (1),  contains  terms  which  result  from  an  analytic  determination  of 
the  ci rcumferenti al  flux  vector,  given  the  assumption  of  axi symmetric  flow  and 

constant  angular  velocity.2  The  viscous  terms  are  contained  in  the  vector  S 
which  is  seen  to  have  variation  in  the  c  direction  only.  This  is  representa¬ 
tive  of  the  thin-layer  approximati on. 

2.  SOLUTION  TECHNIQUE 

Equation  (1)  is  solved  using  the  Bean  and  Warming3  implicit  approximately 
factored  finite  difference  scheme  which  uses  central  differencing  in  both  £ 
and  ;  directions.  Code  improvements  have  been  made  to  include  a  variable  time 
step,  numerical  smoothing  based  on  local  solution  gradients  and  code  vectori- 
zation.1'  The  Beam-Warming  implicit  algorithm  has  been  used  in  various 
applications3'6  for  the  equations  in  general  curvilinear  coordinates.  The 
algorithm  is  first-order  accurate  in  time  and  second-  or  fourth-order  accurate 
in  space.  The  equations  are  factored  (spatially  split),  which  reduces  the 
solution  process  to  one-dimensional  problems  at  a  given  time  level.  Central 
difference  operators  are  employed  and  the  algorithm  produces  block  tridiagonal 
systems  for  each  space  coordinate.  The  main  computational  work  is  contained 
in  the  solution  of  these  block  tridiagonal  systems  of  equations. 

To  suppress  high  frequency  components  that  appear  in  regions  containing 
severe  pressure  gradients  e.g. ,  shocks  or  stagnation  points,  artificial  dissi¬ 
pation  terms  are  added.  In  the  present  application,  a  switching  dissipation 
■odel  is  used  which  is  a  blend  of  second-  and  fourth-order  dissipation  terms. 
This  model  uses  a  fourth-order  dissipation  in  smooth  regions  and  switches  to  a 
second-order  dissipation  in  regions  containing  high  pressure  or  density 
gradients.  Incorporation  of  this  dissipation  model  has  resulted  in  an 
ioprovement  in  the  quality  of  the  results  and  has  made  the  code  more  robust. 
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Tho  developed  ax  1  symmetric  code  uses  a  unique  flow  field  segmentation 
procedure  to  compute  the  full  flow  field  over  a  projectile  or  a  missile 
including  the  base  region.  This  procedure  preserves  the  sharp  corner  at  the 
base.  The  details  of  this  development  can  be  found  in  Reference  6.  For  the 
computation  of  turbulent  flows,  the  two-layer  algebraic  Baldwi n-Lomax  turbu¬ 
lence  model  7  is  used  over  the  body.  A  second  algebraic  turbulence  model  which 
is  based  on  a  simple  exchange-coefficient  concept  is  used  in  the  base  region.8 


Ill.  MODEL  GEOMETRY  AND  COMPUTATIONAL  GRID 

The  external  configuration  of  the  M3  2  5 ,  excluding  the  base,  is  similar  to 
the  M403A1  shown  in  Figure  1.  The  features  of  this  projectile  which  have  not 
been  modeled  exactly  are  the  meplat  on  the  fuze  and  the  rotating  band  near  the 
base.  The  rotating  band  was  eliminated  and  the  meplat  was  modeled  as  a 
hemisphere  cap.  The  computational  model  is  shown  in  Figure  2  and  consists  of 
a  2.34  caliber  nose,  a  2.7  caliber  cylindrical  section,  and  a  0.26  caliber  8° 
boattail.  The  ogive  contours  were  matched  as  well  as  the  undercut  on  the 
cylindrical  section. 

The  current  problem  of  interest  is  the  effect  of  the  different  base 
geometries  on  the  overall  projectile  aerodynamics.  Figure  3a  and  3b  show  the 
standard  and  dome  base  configurations  respectively.  The  standard  base  is  a 
combination  of  aluminum  and  steel  and  contains  a  base  cavity  which  is  charac¬ 
terized  as  a  flat  surface.  The  PIP  configuration  is  an  all  steel  base  and  is 
characteri zed  as  a  domed  surface.  The  cavity  volume  is  also  significantly 
larger  for  the  dome  configuration. 

The  solution  technique  requires  the  discretization  of  the  entire  flow 
region  of  interest  into  a  suitable  computational  grid.  The  grid  outer  bound¬ 
ary  has  been  placed  at  2.3  body  lengths  upstream  and  surrounding  the  projec¬ 
tile.  The  downstream  boundary  was  placed  at  2  body  lengths.  Since  tne 
calculations  are  in  the  subsonic/transonic  regime  the  computational  boundaries 
must  extend  out  beyond  the  influence  of  the  body.  This  ensures  that  the 
boundary  conditions  specified  in  the  flow  code  are  satisfied. 

Figure  4  shows  a  grid  generated  for  the  standard  configuration.  The  grid 
consists  of  225  points  in  the  streanwise  direction  and  50  points  in  the  normal 
direction.  This  is  broken  down  into  two  sections:  a  body  region  and  a  base 
region.  The  surface  points  for  each  region  are  selected  using  an  interactive 
design  program.  Each  grid  section  is  then  computed  separately  using  a  hyper¬ 
bolic  grid  generation  program.9  There  are  165  points  along  the  projectile 
surface,  including  60  points  along  the  afterbody.  The  normal  distribution  of 
points  in  base  region  consists  of  50  points  along  the  base  cavity.  Longitu¬ 
dinally,  the  base  region  is  matched  point  for  point  with  the  60  surface  points 
on  the  afterbody.  An  expanded  view  of  the  base  grid  is  shown  in  Figure  5a. 
The  generally  flat  sections  on  the  standard  base  enabled  a  grid  to  be  routine¬ 
ly  generated.  However,  due  to  the  extreme  concavity  the  grid  for  the  dome 
base  (Figure  5b)  required  an  increase  in  the  smoothing  values  used  by  the 
hyperbolic  grid  generator,  as  well  as  the  addition  of  a  grid  averaging  techni¬ 
que. 
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IV.  RESULTS 


Numerical  computations  have  been  made  for  both  t.he  standard  and  the  dome 
base  configurations  for  the  range  of  Mach  numbers  from  M  =  0.80  to  1.8  and  at 
zero  angle  of  attack.  Solutions  were  marched  in  time  until  the  steady  state 
results  were  achieved.  Atmospheric  flight  conditions  were  used. 

A  few  qualitative  results  are  presented  next.  Figures  6  and  7  show  the 
Mach  number  contours  for  the  dome  base  configuration  for  M  =  0.98  and  1.1, 
respectively.  Figure  6  shows  the  flow  expansion  at  the  ogive  corner,  a  shock- 
wave  on  the  cylinder,  and  an  expansion  at  the  boattail  corner.  In  addition,  a 
shock  system  exists  downstream  of  the  base  corner.  As  the  Mach  number  is 
increased  to  1.1,  the  flow  pattern  has  changed.  For  this  Mach  number,  Figure 
7  shows  a  bow  shock  system  in  front  of  the  nose  of  the  projectile  as  well  as 
the  expansions  at  the  ogive  and  boattail  corners.  The  expansion  at  the  base 
corner  is  much  more  pronounced  and  is  followed  by  the  recompression  shock 
downstream  of  the  base. 

Figure  8  shows  that  the  particle  traces  at  Mach  number,  M  =  1.1  and  a  =  0. 
It  shows  the  flow  expanding  at  ogive  and  boattail  corners.  In  addition,  the 
reci  rcul  atory  flow  in  the  base  region  is  evident.  The  recirculation  region 
extends  to  about  one  caliber  downstream  of  the  base  corner.  Expanded  views  of 
the  particle  trace  plot  in  the  base  region  are  shown  in  Figures  9a  and  9b, 
respectively,  for  the  standard  and  dome  base  configurations.  As  shown  in 
Figure  9a  for  the  standard  base,  the  back  flow,  upon  reaching  the  cavity 
follows  the  contour  of  the  cavity  and  leaves  the  cavity  pushing  the  flow 
upwards.  This,  in  turn,  creates  the  secondary  bubles  seen  near  the  base 
corner.  The  shear  layer  leaving  the  base  corner  is  displaced  upwards  weaken¬ 
ing  the  expansion  at  the  base.  It  clearly  shows  the  primary  bubble  and  the 
secondary  bubbles.  Figure  9b  shows  the  particle  traces  for  the  dome  configur¬ 
ation.  The  flow  again  follows  the  contour  of  the  cavity  and  upon  leaving  the 
dome  cavity,  is  almost  parallel  to  the  streamwise  direction.  This  flow,  thus, 
has  less  effect  on  the  free  shear  layer  and  doesn't  weaken  the  expansion  at 
the  base  corner  as  much  compared  to  the  standard  base.  The  net  effect  is  that 
the  size  of  the  primary  bubble  for  the  dome  base  is  slightly  smaller  than  that 
for  the  standard  base.  The  reattachment  point  is  therefore  closer  to  the  base 
and  should  result  in  lower  base  pressure  or  higher  base  drag. 

Since  the  entire  f 1 ow field  over  the  projectile  including  the  base  region 
is  computed,  all  three  drag  components  can  be  determined  and  thus,  the  total 
aerodynamic  drag  can  be  obtained.  As  stated  earlier,  solutions  were  marched 
in  time  until  the  steady  state  results  were  obtained.  One  of  the  ways  to 
.check  for  convergence  is  to  look  at  the  time  history  of  the  aerodynamic  drag. 
Shown  in  Figure  10  is  the  time  history  for  all  the  drag  components  as  well  as 
the  total  drag  at  M  =  1.1  for  the  standard  base.  As  seen  here,  the  solution 
has  converged  in  2000  time  iterations.  The  viscous  drag  and  pressure  drag 
components  have  converged  in  less  than  500  iterations.  However,  as  expected, 
it  takes  a  lot  longer  for  the  base  region  flow  to  converge.  The  base  flow  is 
oscillating  in  the  transient  stage  but  the  final  converged  result  indicate 
that  the  flow  is  steady.  This  result  is  typical  of  the  high  transonic  Mach 
number  (M  >  1.1)  runs  for  both  base  configurations.  At  low  transonic  Mach 
numbers  (M  <  1.1),  some  unsteadiness  in  the  flowfield  has  been  observed.  A 
i  typical  plot  for  such  a  case  is  shown  in  Figure  11  for  Mach  =  0.8.  As  seen  in 

\  this  figure,  the  base  drag  as  well  as  the  pressure  drag  both  indicate  a 
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periodic  unsteady  f  1  ov/  behavior.  The  drag  values  for  these  cases  have  been 
taken  as  the  averages  of  the  amplitude  extremes. 

Figure  12  shows  the  pressure  drag  as  a  function  of  Mach  number.  The  pres¬ 
sure  drag  rise  in  the  transonic  speed  regime  is  clearly  predicted.  The  dif¬ 
ference  in  pressure  drag  between  the  two  projectiles  (standard  base  and  dome 
base)  is  very  small  at  high  transonic  Mach  numbers  and  gets  larger  (=  15%)  at. 
the  low  transonic  Mach  numbers  (M  =  0.9).  The  viscous  drag  component  is  shown 
in  Figure  13.  This  component  of  the  total  drag  is  rather  small  and  the 
difference  in  viscous  drag  due  to  the  base  configurations  is  negligible.  The 
difference  in  the  shapes  of  the  cavities  is  expected  to  impact  the  base  drag 
which  is  shown  in  Figure  14.  The  difference  in  base  drag  between  the  standard 
base  and  the  dome  base  is  larger  at  the  high  transonic  Mach  numbers  (1.1  <  M  < 
1.1)  as  well  as  at  low  transonic  speeds  (M  <  0.94).  The  difference  is  very 
small  near  M  =  0.97.  In  addition,  the  base  drag  is  higher  for  the  dome  base 
than  the  standard  base  at  high  transonic  Mach  numbers  (0.97  <  M  <  1.1)  while 
the  reverse  is  true  at  low  transonic  Mach  numbers  (M  <  0.97).  The  expected 
rise  in  base  drag  at.  transonic  speeds  is  clearly  predicted.  Comparison  of  the 
total  aerodynamic  drag  is  shown  in  Figure  19.  As  shown  in  this  figure,  the 
difference  in  drag  is  very  small  near  M  =  0.97  and  is  somewhat  larger  at.  high 
transonic  speeds  (1.1  <  M  <  1.9)  as  well  as  at  low  transonic  speeds  (M  < 
0.92).  This  plot  also  shows  the  range  data  for  both  base  configurations.  The 
overall  comparison  of  the  computed  drag  data  with  the  range  data  is  fair.  As 
seen  in  the  range  data,  the  dome  base  has  higher  drag  especially  at  higher 
transonic  Mach  numbers  and  this  trend  is  seen  in  the  computed  results  also. 


V.  CONCLUDING  REMARKS 

Numerical  computations  have  been  made  for  a  projectile  with  two  base 
cavity  configurations  at  transonic  speeds.  Computed  results  show  differences 
in  the  qualitative  features  of  the  base  reqion  flowfield  between  the  two  base 
cavities.  The  drag  components  and  the  total  drag  have  been  computed.  Changes 
in  the  base  cavity  configuration  have  been  found  to  affect  the  aerodynamic 
drag.  Oifferences  in  drag  of  between  0  to  19%  have  been  predicted.  The  dome 
hase  configuration  produces  a  higher  drag  especially  at  higher  transonic 
speeds  (1.1  <  M  <  1.9)  and  qualitatively  agrees  with  the  trend  of  the  experi¬ 
mental  range  data.  At  lower  transonic  spe.ds  (M  <  0.97),  the  reverse  is  true. 
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Figure  5b.  Base  region  grid  for  the  dome  base. 
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LIST  OF  SYMBOLS 


a  =  speed  of  sound 

D  =  body  diameter 

e  =  total  energy  per  unit  vol ume/p^a^ 

£,  G,  =  flux  vector  of  transformed  Navier-Stokes  equations 
h  =  time  step 

H  =  n-invariant  source  vector 

J  =  Jacobian  of  transformation 

i'  =  Mach  number 

p  =  press  ure/pjs* 

q  -  dependent  variables 

P.e  =  Reynolds  number, 

S  =  viscous  flux  vector 

t  =  physical  time 

u,  v,  w  =  Cartesian  velocity  component s/a^ 

x,  y,  i  -  physical  Cartesian  coordinates 

Greek  Symbols 

a  =  angle  of  attack 

u  =  coefficient  of  vi scosity/u^ 

s ,  n ,  ;  =  transformed  coordinates  in  axial,  ci rcumf erenti al  and  radial 

di rections 

p  =  density/p^ 
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